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Abstract 

This paper derives the total power or energy loss rate generated in the form of gravitational waves 
by an inspiralling compact binary system to the five halves post-Newtonian (2.5PN) approximation 
of general relativity. Extending a recently developed gravitational-wave generation formalism valid 
for arbitrary (slowly-moving) systems, we compute the mass multipole moments of the system 
and the relevant tails present in the wave zone to 2.5PN order. In the case of two point-masses 
moving on a quasi-circular orbit, we find that the 2.5PN contribution in the energy loss rate is 
entirely due to tails. Relying on an energy balance argument we derive the laws of variation 
of the instantaneous frequency and phase of the binary. The 2.5PN order in the accumulated 
phase is significantly large, being grossly of the same order of magnitude as the previous 2PN 
order, but opposite in sign. However finite mass effects at 2.5PN order are small. The results of 
this paper should be useful when analyzing the data from inspiralling compact binaries in future 
gravitational-wave detectors like VIRGO and LIGO. 
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I. INTRODUCTION 



Compact binaries in their late stage of evolution are very relativistic systems in which the 
two compact objects (neutron stars or black holes) orbit around eachother with velocities 
as large as 30 % of that of light. The gravitational radiation these systems emit during 
the inspiral phase preceding the coalescence of the two objects is expected to be routinely 
analyzed in future detectors like LIGO and VIRGO (see for reviews). Hundreds to 

tens of thousands of gravitational- wave cycles should be monitored in the detectors' sensi- 
tive frequency bandwidth. The combination of high orbital velocities and a large number 
of observed rotations, together with the fact that the emitted waves are highly predictable, 
implies that high-order relativistic (or post-Newtonian) effects should show up in the grav- 
itational signals observed by VIRGO and LIGO fl fl S Q, B S Q- Alternatively, this 
means that high-order post-Newtonian effects should be known in advance so that they can 
be included in the construction of theoretical filters (templates) to be cross correlated with 
the outputs of the detectors. 

The relevant model for describing most of the observed inspiral phase is a model of two 
point masses moving on a circular orbit. Radiation reaction forces tend to circularize very 
rapidly the orbit. On the other hand, point masses can be used in the case of non-rotating 
and (initially) spherically symmetric compact objects up to a very high precision IjJ. This is 
due to a property owned by general relativity of "effacing" of the internal structure. Even in 
the case of stars with intrinsic rotations the dynamics of the binary is likely to be dominated 
by post-Newtonian gravitational effects jl^ . 

High-order post-Newtonian effects that are measurable are mainly those affecting the 
orbital phase evolution of the binary, which in turn is determined using a standard energy 
balance argument by the total power emitted in the form of gravitational waves by the 
system at infinity, or total luminosity in the waves. To what level in a post-Newtonian 
expansion should we know the gravitational luminosity (or energy loss) in order to guar- 
antee an optimal detection of the signal (given some power spectral density of the noise in 
a detector) is still unclear, but the theory of black-hole perturbations can be used to gain 
insights in this problem. Black-hole perturbations, which deal with the special case of a 
test mass orbi ting a massive black hole, have been recently the focus of intense activity 



0Q,Q,0,lii 



It emerges from this field that neglecting even a so high approximation 
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as the third post-Newtonian (3PN) one, i.e. neglecting the relativistic corrections in the 
luminosity which are of relative order (or below)j_is likely to yield unacceptable system- 



li 



181 19^ . This shows how relativistic are 



atic errors in the data analysis of binary signals 
inspiralling compact binaries, as compared for instance to the binary pulsar for which the 
Newtonian approximation in the luminosity (Einstein quadrupole formula) is adequate. The 
post-Newtonian theory is presently completed through the second post-Newtonian (2PN) 
approximation, i.e., through relative order (both in the waveform and in the associated 
ener gy l oss ) . Two computations were performed to this order, one by Blanchet, Damour and 
yer |2n[ I21I based on a post-Minkowski matching formalism, and one by Will and Wiseman 
22 1 using an approach initiated by Epstein and Wagoner j23[ and generalized by Thorne 
The common result of these two computations for the energy loss was summarized in 
jl^ . and the waveform can be found in 12^" 



Ref. 



In the present paper we develop the post-Minkowski matching formalism one step beyond 



the work of Refs. [2ij|,|2l| by computing the 2.5PN order in the energy loss of an inspiralling 
binary. This entails extending both Ref. 3] on the general formalism valid for an arbitrary 



(slowly-moving) source, and Ref. j2l| dealing with the specific application to the binary. 
The computation of the waveform of the binary to 2.5PN order (the square of which should 
give back the 2.5PN energy loss) will be left for future work. 

The post-Minkowski matching formalism is a wave generation formalism which is espe- 
cially suited for "semi relativistic" sources whose internal velocities can reach 0.3c at most 
(say), as in the case of inspiralling compact binaries (see [2^] for a review). The formalism 
combines (i) an analytic post-Minkowskian approximation scheme for the computation of 
the gravitational field in the exterior of the source where multipole expansions can be used to 
simplify the problem, (ii) a direct post-Newtonian approximation scheme for the resolution 
of the field equations inside the near-zone of the source, and (iii) an asymptotic matching 
between both types of solutions which is performed in the exterior part of the near zone. The 
necessicity of using a post-Minkowskian approximation scheme first is because its validity 
extends up to the regions far away from the system where the observer is located, contrar- 
ily to the post-Newtonian approximation whose validity is limited to the near zone. The 
exterior field is computed using an algorithm developed in Ref. which set on a general 
footing previous investigations by Thorne i3] and Bonnor The implementation of the 
wave generation formalism (i)-(iii) was done at first with IPN accuracy in Ref. [2^ which 
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obtained the IPN correction terms in the mass-type quadrupole moment of the source (and 
in fact in all the mass-type multipole moments). The dominant nonlinear contribution in the 
radiation field was added in Ref. js^ and shown to be due to the contribution in the wave 
zone of the well-known "tail" effect. The inclusion of this nonlinear contribution pushed the 
accuracy of the formalism to 1.5PN order in the energy loss. The 2PN precision both in the 
energy loss and waveform was reached in Ref. Q| where the second (2PN) correction terms 
in the mass-type multipole moments and the first (IPN) ones in the current-type multipole 
moments were obtained. An equivalent expression of the IPN current-type moments had 
been derived earlier , however in a different form which is not amenable to generalization 
to higher order (but turned out to be useful in applications |2ll|). 

The main result of the present paper is the expression of the 2.5PN-accurate energy loss 
by gravitational radiation from a general (semi relativistic) source, and from an inspiralling 
compact binary. In the latter case of application, the 2.5PN contribution in the energy loss 
is found to be entirely due to tails in the wave zone (this is like the 1.5PN contribution) , 



00,13. 



and to reduce in the test-body limit to the known result of perturbation theory 

With the energy loss one can derive the laws of variation of the inspiralling binary's 
orbital frequency and phase using an energy balance equation. However note that this is a 
weak point of the analysis because the latter energy balance equation has been proved to 



hold only at the Newtonian order (see [Sj, |33|, |3J, |35|, 36 1 for general systems, and ^37] for 
binary systems), and more recently at the IPN order 2S) It is also known to hold for 
the specific effects of tails at 1.5PN order To prove this equation at the 2.5PN order 

as we would need below, one should in principle obtain the equations of motion of the binary 
up to the very high 5PN order (or order c^^°) beyond the Newtonian acceleration. Indeed 
the radiation reaction forces which are responsible for the decrease of the binding energy of 
the binary are themselves dominantly of order 2.5PN beyond the Newtonian term. The only 
method which is available presently in order to deal with this problem is to assume that the 
2.5PN-accurate radiation reaction forces are such that there is exact agreement between the 
loss of 2.5PN-accurate binding energy of the binary (as computed from the D amour- Deruelle 



equations of motion 



37| ) and the 2.5PN-accurate energy flux we shall compute below. This 



assumption is verified at the 1.5PN order and sounds reasonable, but will have to be justified 
in future work. 

As an indication of the quantitative importance of the 2.5PN approximation in the or- 
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bital phase of the binary, we compute the contribution of the 2.5PN term to the number 
of gravitational-wave cycles between the entry and exit frequencies of some detector. Es- 
sentially we find that the 2.5PN approximation is of the same order of magnitude as the 
2PN term (computed in but opposite in sign. In the case of two neutron stars of mass 

IAMq and of the frequency bandwidth [10 Hz, 1000 Hz], the 2PN term contributed to +9 



units to the total number of cycles [1^. We find here that the 2.5PN term contributes to 
— 11 cycles in the same conditions. This shows the importance of the 2.5PN approximation 
for constructing accurate theoretical templates. [However, we find that the contribution of 
the finite mass effects in the 2.5PN term (which cannot be obtained in perturbation theory) 
is numerically small.] 

In the present paper we shall make a thorough investigation (see Sections II, HI and 
IV below) of all the relativistic corrections in the multipole moments of the system which 
contribute to the 2.5PN-accurate energy loss. However when we are interested only in the 
application to inspiralling compact binaries, this investigation can be seen a posteriori to 
be unnecessary. Indeed the orbit of an inspiralling binary is circular, and we shall prove 
in this case that the 2.5PN relativistic corrections in the multipole moments give in fact 
no contribution in the energy loss. As we said above, the only contribution is that of the 
tails present at this order (see Section VI). A simple argument (concerning the result at 
2.5PN order of a contracted product of tensors made of the relative separation and velocity 
of the bodies) could be used beforehand to see that this is true. But because inspiralling 
compact binaries may not constitute the only sources for which the 2.5PN approximation 
in the energy loss is needed, or simply because one may need in the future to consider the 
case of a binary moving on an eccentric orbit, we have chosen in this paper to compute 
systematically all the terms which enter the 2.5PN-accurate energy loss for general systems. 
This permits us to show explicitly that all the terms but the tail terms give zero in the 
energy loss for (circular) inspiralling binaries. The simple argument mentioned above may 
be used in future work to simplify the investigation of higher post-Newtonian orders. 

The plan of this paper is as follows. In the following section II and in the next one HI we 
follow step by step the derivation done in Ref. 2^ of the near-zone gravitational field and 
the corresponding matching equation, and show how this derivation can be extended to the 
2.5PN order. In section IV we obtain the explicit 2.5PN corrections arising in the mass-type 
multipole moments. Section V deals with the derivation of the energy loss formula valid 
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for general systems (however some coefficients are left un-specified in the formula). Finally 
these results are applied in Section VI to inspiralling compact binaries. Appendix A derives 
a useful integration formula and Appendix B presents a relevant summary of the 2.5PN 
equations of motion. 



Throughout this paper we refer to Ref. as paper I and to Ref. 



2l| as paper II. 



II. THE GRAVITATIONAL FIELD IN THE NEAR ZONE 

Following the plan of paper I we first investigate the gravitational field generated by a 
slowly moving isolated source in its near zone, which is defined in the usual way as being 
a zone whose size is of small extent with respect to a typical wavelength of the emitted 
radiation. Two distinctive methods are used. The first one is a direct post-Newtonian 
iteration (speed of light c +00) of the field equations inside the source, and is valid all- 
over the near zone. The second method consists of re-expanding when c +00 a solution 
of the vacuum field equations obtained by means of the multipolar and post-Minkowskian 



iteration scheme of Ref. j27|, and is valid only in the exterior part of the near zone. 

We denote the small post-Newtonian parameter by e ~ f/c, where f is a typical velocity 
in the source (e.g., the relative orbital velocity of the two bodies in the case of a binary 
system). A remainder term of order 0(6:") is abbreviated by 0{n). In a vector or a 
tensor B^"^ , the remainder term is denoted by 0{n,p) or 0{n,p,q), by which we mean a 
term of order 0{n) in or of order 0{p) in A'' or B^^ or and of order 0{q) in 
B^^ . [Greek indices range from to 3, and Latin indices range from 1 to 3.] Most of the 
notations used here are like in paper I. 



A. The inner gravitational field 

T"'^ are the contravariant components of the stress-energy tensor (with dimension of an 
energy density) of the material source in some inner coordinate system (x, t). The densities 
of mass a, of current and of stress aij in the source are defined by 

T^OO _|_ n^ii 

o = t ' (2-la) 

a^ = — , (2.1b) 

c 
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Oij = , (2.1c) 

where T** denotes the spatial trace T^SijT^^ . These definitions are such that a, ai, Uij have 

a finite non-zero hmit as c — > +00. Prom these matter densities one defines the retarded 
potentials 

V = -47rG'n;^V , (2.2a) 

Vi = -47rG'n^Vi , (2.2b) 

1 



(2.2c) 



where G is Newton's constant, and where denotes the retarded integral operator 

(□hVKx,^)^-^/// ^^/(x',t-l|x-x'|). (2.3) 

Contrarily to the sources of V and Vi which are of compact support: UV — —AnGa and 
nVi — —AnGai where □ is the d'Alembertian operator, the source of the potential Wij is not 
of compact support: OiWij = —AirGaij — diVdjV + ^5ijdkVdkV. Indeed we have included in 
Wij the stress density of the (Newtonian) gravitational field itself since it is of the same order 
as (Tjj when c ^ +00. [Note that paper II used the notation Wij for a closely related but 
different potential; here we do not follow paper II but stick to the notation of paper I.] To 
Newtonian order the densities and potentials so defined satisfy the equations of continuity 
and of motion 

dta + diai = 0(2) , (2.4a) 
dtai + djGij = adiV + 0(2) . (2.4b) 

Prom these dynamical equations one deduces the differential identities 

dtV + diVi = 0(2) , (2.5a) 
dtVi + djWij = 0(2) . (2.5b) 

With the introduction in paper I of the retarded potentials V , Vi, Wij, a simple expression 
of the gravitational field h""^ inside the source which is valid to some intermediate accuracy 
0(6, 5, 6) was written, namely 



/i™ = -^V + ^{Wii - 2V^) + 0(6) , (2.6a) 



h'' = --^V, + 0{5) , (2.6b) 



c 



h'^ = -^Wi, + 0(6) . (2.6c) 



The field variable is h"^ = \/—gg°'^ — rj"'^, where g and g"'^ are the determinant and inverse 
of the usual covariant metric gajj, and where 77"^ is the Minkowski metric (signature — 1-++). 
Note the important fact that there are no explicit terms in Eqs. ()2.6|) involving powers of 
which are "odd" in the post-Newtonian sense (e.g., a term of order ~ in or 
h^^). This is because we have kept the potentials V, Vi and Wij in retarded form, without 
expanding when c +00 the retardation they contain. The "odd" terms in Eqs. (j2.6|) could 
be easily computed using the post-Newtonian expansions of the retarded potentials as given 
by Eqs. (jOj) below. 

Paper I iterated the inner field ()2.6p from this intermediate post-Newtonian order to the 
next order with the result that the field to the higher precision 0(8, 7, 8) could be written 

as 

+ 0(8,7,8) (2.7) 



where A and A denote some explicit combinations of the inner gravitational potentials V , 
Vi, Wij and their derivatives. First we have 

\{y, W) = l + \v- -AWu - V^) (2.8) 

which represents in fact the post-Newtonian expansion of (minus) the determinant of the 
metric, with terms 0(6) suppressed. Second the components of A read 



-rOO, 
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K'-'iy.W) = --dkVdkV + -^-VdtV - 2VudtdkV 

-Wkmdl^V + lidtVy + IdkVUdkVm + 3dM 

o Z 

. 20.ve.w^^ - Ivo^vo^v] . ,2.0a, 

A^'iV, W) = ^ l^dkVid.Vk - dkV^ + ^dtVd.v'^ , (2.9b) 

3. 



-dkVAv, + 2d^iyAVj) - -6,,{dtvy - s^.d^vdtVk 
1 

2' 



+]:5ijdkVm{dkVm - dmVk) \ , (2.9c) 



and represent the expansion of the effective nonhnear gravitational source of Einstein's 
equations in harmonic coordinates, with 0(8, 7, 8) terms suppressed. The overbar on A and 
A reminds us that these quantities are only determined up to a certain post-Newtonian 
order. The (approximate) harmonic coordinate condition is 



As we shall see the post-Newtonian accuracy of the inner field fl2.7|l - ()2.1()|l is sufficient for 
our purpose. 

B. The external gravitational field 

In the exterior we use a solution of the vacuum field equations which has in principle 
sufficient generality for dealing with an arbitrary source of gravitational radiation. This 
solution is given as a (nonlinear) functional of two infinite sets of time-varying multipole 
moments, M^lt) and 5*2, (t). The index L carried by these moments represents a multi-index 
formed with i spatial indices: L = 11X2 . . .ie. The "order of multipolarity" i goes from zero 
to infinity for the "mass-type" moments ML{t) and from one to infinity for the "current- 
type" moments Siit). The moments Ml and Sl are symmetric and tracefree (STF) in their 
i indices. The mass monopole M is simply the total mass of the source or ADM mass, the 
mass dipole Mj is the position of the center of mass (in units of total mass), and the current 
dipole Si is the total angular momentum. M, Mj and Si are constant. Furthermore we shall 
choose Mi = by translating the spatial origin of the coordinates to the center of mass. 

Some external potentials playing a role analogous to the inner potentials but differing from 
them both in structural form and numerical values are introduced. First the potentials V^'^^*, 



yext g^j^^ Y^xt g^j^g given by their multipolar series parametrized by the multipole moments 
Miit) and Siit) : 




(2.10) 




(2.11a) 





(2.11b) 
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ext 



t=2 

oo 



-d 



L-2 



2e 



e=2 



i\ £ + 1 



da 



L-2 



r 

t - - 

c 



(2.11c) 



The notation 9l is shorthand for a product of partial derivatives, (9l = di^di^ ■ ■ - dig where 
di = d/dx\ In a similar way = . . . = dad^-i, and so on. The superscript 

(n) indicates n time derivatives, and the indices in brackets are symmetrized. The potentials 
fj2.11|) satisfy the source-free d'Alembertian equation and the (exact) differential identities 



ext 








(2.12a) 
(2.12b) 



Furthermore V^Y^ is tracefree: V^f^^ = 0. Having defined these potentials one introduces a 



more complicated potential W^j^^ by the formula 



ext 



l/.ft + FP 



B=0 



□ 



R 



|^_5^V^ext^.^ext ^ l^.^.^^^ext^^^ext^ 



(2.13) 



The second term appearing here involves the retarded integral (j2.3j) but regularized by 



means of the analytic continuation process defined in Ref. [27[. One multiplies the source 
by where r = |x| and i? is a complex number. Then applying the operator defines 
a function of B which can be analytically continued all-over the complex plane deprived 
of the integers, but admitting there a Laurent expansion with only some simple poles j^]. 
The looked-for solution is equal to the finite part at -B = (in short FPb=o) or constant 
term ~ in the Laurent expansion when — > 0. This regularization process is made 
indispensable by the fact that we are looking for solutions of the wave equation in the form 
of multipole expansions, which are valid only in the exterior of the source and are singular 
when considered formally inside the source. 

With the external potentials V^^^*, Vf^'^ and W^J^^ one has a result similar to Eq. ()2.(ij] . 
namely the expression of the external field h^^!^^ (where the notation "can" means that we 
are considering specifically the "canonical" construction of the external field as defined in 
Sect. 4.3 of |,27,j) up to the intermediate accuracy 0(6,5,6) : 



h 



00 



-4 V' + 4 (w^r* - 2(K'="*)') + 0(6) , 



Oi 



vr' + 0(5) 



(2.14a) 
(2.14b) 
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/^L = -^W:f + 0{Q). (2.14c) 

Because this expression has the same form as Eq. ()2.6p . the nonhnearities in the exterior 
field will have in turn the same form as Eq. ()2.9j) . Up to order 0(8,7,8) we obtain the 
"canonical" field 

= GKlw + G"€Lm + FPb=o [r^T\V^^\ ly-*)] + 0(8, 7, 8) , (2.15) 
satisfying the (exact) harmonic gauge condition 

duKL = . (2.16) 

In the last term of (j2.15p the symbol FP b=o and the regularization factor have the same 
meaning as in ()2.13|) . The effective nonlinear source A^'^ is given by Eqs. ()2.9|) but expressed 
with the external potentials V^'^^*^ Vf^^ and W^J^^ instead of the inner potentials V, Vi and Wij 
appearing in ()2.7|) . The first term Gh^^^^-^^ in Eq. ()2.15p is linear in the external potentials, 

A 

(2.17a) 
(2.17b) 
(2.17c) 

This term is the solution of the linearized (vacuum) equations on which is based the post- 

n 

Minkowskian algorithm for the construction of the canonical metric in Ref. |27|]. The powers 
of in Eqs. fITTTjl are such that h^l^^^-^ = 0(2,3,4). Finally the term G'^q'^^^^^^ in (|TT^ 
is a particular solution of the source-free wave equation which has to be added in order 
that the harmonic gauge condition (j2.16j) be satisfied. It was proved in the Appendix A of 
paper I that this term is of order 

GV:n{2)=0(7,7,7), (2.18) 

and thus can be safely neglected if we are interested in the mass multipole moments to 2PN 
order only. In the present work investigating the 2.5PN order the term ()2.18|) cannot a 
priori be neglected. However we shall see that the 2.5PN order is needed only in the sum 
of the 00 component and the spatial trace of this term, i.e. G'^{q^^^i^2) + ?can{2))- Relying 
on our previous papers we know that this sum is made of some retarded waves of the type 



'^can(l) — 


4 

~^ 


"^ext 


^^can(l) = 


4 

~^ 


^ i 


Gh'^ = 
can(l) 


4 


Tyext 
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dL{r~^X{t — r/c)) with scalar or dipolar multipolarity only = or 1). Furthermore we 
know that the dependence on of such a wave is 0(5 + + £3 ~ ^) where ii and £2 
the multipolarities of the two interacting moments composing the wave (we use the same 
notation as e.g. the Appendix A of paper I; see also after Eq. below). With £ = or 

1 a term of order 0(7) necessarily has £^ + £3 = 2 or 3. The term in question is made of 
the product of the mass M with the quadrupoles Mij and Sij, or of the product of the mass 
dipole Mi with the quadrupole Mij. The first possibility is excluded because My and S^j 
are STF (thus no scalar or dipole wave with no free index can be formed), and the second 
possibility does not exist in a mass centred-frame where Mj = 0. So we have proved that in 
a mass-centred frame we have 

G''(A2)+Cn(2)) =0(8) , (2.19) 

which is all what will be needed in the following. 



III. THE MATCHING BETWEEN THE INNER AND OUTER FIELDS 
A. Relations between inner and outer potentials 

Our matching requirement is that there exists a change of coordinates valid in the external 
near zone and transforming the inner gravitational field /i"^(x) given by Eq. (j2.7|) into the 
outer field /^"^(xcan) given by Eq. ()2.15|) . Let this change of coordinates be 

xi^,^{x) = x^ + V^{x) (3.1) 

where are the inner coordinates used in section II. A and x^f^^ the outer coordinates used 
in section II. B. The vector ip^ has been shown in paper I to be of order 

^'^ = 0(3,4). (3.2) 

Because the two coordinate systems and x^f^^ are harmonic [at least approximately ; see 
Eqs. (^J(l\i and ^U^ ]. and because /i"^ = 0(2,3,4) in addition to Eq. we have 

Dip^" = 0{7,8) . (3.3) 

The matching equations consistent with the order 0(8, 7, 8) read from paper I as 

hcLi^) = h'^'ix) + dip''' + 2/i°'^9^(^° - 9^(/i°V^) + ip' d, ip' + 0(8) , (3.4a) 

12 



h'U^) = h^'{x) + dv'^ + 0{7) 
h^U^) = h'={x) + d^'^+0{8) 



(3.4b) 
(3.4c) 



where both the inner and outer fields are expressed in the same (inner) coordinate system 
x^. We denote by dip'^'^ the hnear part of the coordinate transformation, 



(3.5) 



The nonhnear part of the transformation enters to this order only the 00 component ()3.4ap . 
The matching equations were used in paper I first to obtain the relations between the 
external potentials and the multipole expansions of the corresponding inner potentials. The 
insertion of the intermediate expressions of h^''' and h^^^ valid up to 0(6, 5, 6) and given by 
and fTWf yields 



■^ext 



ext 



MiV) + cdtip'' + OiA) , 
MiVi}-jd,ip^ + 0{2) , 

M{W^,) - j[d,ip^ + d.if' - 5,,(9o</^° + dk^')] + 0(2) . 



(3.6a) 
(3.6b) 

(3.6c) 



The script letter Ai refers to the multipole expansion. Note that the first equation is valid to 
post-Newtonian order [see the remainder 0(4)] while the two others are valid to Newtonian 
order only [remainder 0(2)]. The multipole expansions of V, Vi and Wij are given by 



£=0 

1=0 ^■ 



1=0 



iv- (t - 

r \ c 
r \ c 



-W^. (t - 

r *H c 



(3.7a) 
(3.7b) 



-dMiV)d,MiV) + -5,,dkMiV)dkMiV) 



, (3.7c) 



where the (reducible) multipole moments V^, and W^j are 
V^(t) = J d^yiXL j dz5t{z)a{y.,t + z\y.\/c) , 



vf(t) 



1 

df'^XL / dz 5i{z)oi{yi,t + z\-x\/ c) 



(3.8a) 
(3.8b) 



13 



Oi, + ^ (d,Vd,V - \b,,d^Vd^V^ (x, t + z|x|/c) . (3.8c) 

The notation is for the tracefree projection of the product of I spatial vectors x/, = 
Xi^ - ■ ■ x^. The function 5i{z) is given by 



5,{z) 



^^i^T^(l-^^r ; Jjz6e{z) = l. (3.9) 



Some explanations on the expressions ()3.7|) - ()3.9p are in order. First notice that the expres- 
sions of the multipole expansions of the potentials V and Vi whose sources have a compact 
support are quite standard. They can be found in this form in the Appendix B of Ref. j^] 
(but were derived earlier in an alternative form 4l|). Notably the presence of the function 
Se{z) is due to the time delays of the propagation of the waves with finite velocity c inside 
the source. By contrast the expression of the multipole expansion of the potential Wij whose 
source extends up everywhere in space is more complicated. The second term in Eq. (j3.7cp 
ensures that J^{Wij) satisfies the correct equation [deduced from ()2.2cj) ] outside the source, 
namely DM{Wij) = -diM{V)djM{V) + l6ijdkM{V)dkM{V) where the right side agrees 
numerically with —diVdjV + \5ijdkVdkV outside the source. This second term involves the 
multipole expansion M.(y) of the inner potential and not V itself. This is in conformity 
with the use of the regularized operator FP which is defined only when acting on mul- 
tipole expansions (like in Eq. ()2.13|l ). On the other hand the integrand of ()3.8c|l involves the 
non-compact supported source of Wij where appears the potential V itself, i.e. not in mul- 
tipole expanded form. Thus the integrand of ()3.8c|l is valid everywhere inside and outside 
the source. Very far from the source it diverges because of the presence of the product x^ of 
I spatial vectors, behaving like |x|^ at spatial infinity. The well-definiteness of the integral 
results from the presence of the regularization factor |x|^ and the finite part symbol. We 
notice that no ad hoc prescription is necessary in order to obtain the multipole moments in 
a well-defined form even in the case of a non-compact supported source. This is proved in 
paper I. 
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B. The matching equation 



C3 



The relations ()3.6p linking V'''^^ , V^^"", Wif- to the multipole expansions M{V), M{Vi), 
}A{Wij) serve us with re-expressing the non-linearities in the external metric in terms of the 
potentials belonging to the inner metric. The result of paper I is 

F"(V^""*,W^""*) =F"(A<(l^),X(iy)) + m'^" + 0(8,7,8) , (3.10) 

where the components of the tensor Vt^^ are given by 

= 

^oi _ ^ij ^ Q _ (3.11b) 
The only non-zero component of this tensor is the 00 component which is of order 

fi°° = 0(6). (3.12) 

This is because = 0(3), (^9* = 0(4). Both sides of Eq. ()3.10p are multiplied by r^, and 
one applies afterwards the retarded integral and take the finite part at 5 = 0. We find 

FPB=onR^ [r^A'"'(F""*, ly^"*)] = FPB=on:R^ [r^T\M {y),M {W))] + fi'^^ + X^'^ + 0(8, 7, 8) 

(3.13) 

where one must be careful that an extra term X^'^ with respect to paper I arises whose 
components are given by 

= FPB=on^'[^''™°°] -^"^ , (3-14a) 
^oi ^ -^ij ^ Q _ (3.14b) 

The only non-zero component of this term is the 00 component which satisfies 

= ; = 0(7) . (3.15) 

The fact that is a (retarded) solution of the source-free d'Alembert equation is clear 
from its definition (j3.14ap and the main property of the operator FFD]^^ which is the inverse 
of □ (when acting on multipolar sources). The fact that is of order 0(7), that is exactly 
one order in smaller than the order of the corresponding Q^^, is not immediately obvious 
but will be proved below. The quantity was rightly neglected in paper I but has to be 
considered here since it will contribute to the mass moments at the 2.5PN order. 
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We are now in a position to write down a matching equation valid up to the neglect of 
0(8, 7, 8) terms. The equation is obtained by insertion into Eqs. ()3.4|) of both the inner field 
h^'-' given by (|2.7|) and the outer field /i^f^^ given by (j2.15|) . Use is made of the link derived 
in (j3.13p between the external and internal non-linearities. We find that Q^'^ cancels to the 
required order the nonlinear part of the coordinate transformation so that only remains the 
linear part dip'^'^ given by Eq. ()3.5|) . The resulting matching equation (extending the less 
accurate equation (3.35) of paper I) reads as 



can(l) ^can 2) H 



can(l) "can(2) 



C 



A 



-X^^ + a(/9^^ + 0(8,7,8) . (3.16) 

There are two new terms with respect to (3.35) in paper I : G^g^^^-g-) in the left side and 
—X^^ in the right side which both are terms of 2.5PN order. Now, by exactly the same 
reasoning as in paper I one can transform the difference between the two retarded integrals 
in the right side into an explicit multipole expansion parametrized by some moments {t\ 
namely 



1=0 



r 



X^"" + d^'"' + 0(8,7,8) . (3.17) 



These moments are given by 



Tl\t) = FPb=o j d'^^l^l^XL j dz5^{z)T^"'{^,t + z\^\/c) (3.18) 

where T'^" denotes the total stress-energy tensor of the material and gravitational fields (valid 
up to the considered precision) 

rf'-^V, W) = XiV, W)T^' + -^T\V, W) , (3.19) 

167rG 

which is conserved in the sense 

d^T^"^ = 0{?>,A) . (3.20) 

It is important to note that the effective stress-energy tensor t^^ is a functional of the poten- 
tials V , Vi and Wij valid everywhere inside and outside the source. There is no contribution 
of the multipole expansions of the potentials in the final result (see paper I). 
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The left side of the matching equation ()3.17|) is a functional of the original multipole 
moments Ml, Sl parametrizing the exterior metric. On the other hand the right side is 
a functional of the actual densities of mass, current and stress of the material fields in the 
source. To find the explicit expressions of and Sl in terms of these source densities we 
decompose the right side into irreducible multipole moments. Inspection of the reasoning 
done in section IV.A of paper I shows that this reasoning is still valid in the present more 
accurate case. As a result we find 



G'Cn(i) [Ml. Sl] + G'qZ^^,^ = Gh^^^^ [h, Jl] - + + 0(8, 7, 8) 



(3.21) 



where the last term is a linear coordinate transformation associated with the vector = 
^p^^ _|_ ^JJ^^ where uj^ is the same as in Eq. (4.4) of paper I. The linearized metric Gh^^^-^^ in 
the right side takes the same expression as in the left side but is parametrized instead of 
Ml and 5*^ by the (STF) source multipole moments II and Jl given by 



hit) = FPb=o / c^'x|x 



dz 



6i{z)xlT. 



4(2£+ 1) 



c2(£+l)(2£ + 3) 



6i+i{z)xiLdtT.i 



2(2£+ 1) 



5i+2{.z)XijLdlTjij 



(x,t + z|x|/c) , (3.22a) 



JL{t) = FPb=o / d^^\^f j dz 

21+1 



c2(£ + 2)(2£ + 3) 



{z)eab<iiXL~l>acdt^hc 



We have posed 



7^0 _j_ jii 



Eij — T 



C 



{x,t + z\x\/c) (3.22b) 

(3.23a) 

(3.23b) 
(3.23c) 



The equation ()3.21|1 can be solved uniquely for the multipole moments Ml and Sl- To do 
so it suffices to notice that for any gauge term dC,^'^ = d^S,'^ + d'^E,^ — ri^^d\^^ the identity 

- + di''^] + doidM""^ + + dl[dc^ + -5'^(ae°° - d^''^)] = o (3.24) 

2 ■' 2 

holds. This is nothing but the vanishing of the OiOj component of the linearized Riemann 
tensor when computed with g^^^ = d^j,^^, + 8,^^^. Applying ()3.24|) to the gauge term of 
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()3.21|) and using the form of the multipole moment decomposition of ^can(i) finds 



Ml = Il + SIl + 0(6) 
Sl = Jl + 0{4) , 



(3.25a) 
(3.25b) 



where SIl in the mass moments is of 2.5PN order and comes from the decomposition of 
X"" (which we recall is a retarded solution of the wave equation) into multipole moments 
according to 



X 



00 



H 



-dr 



-5h{t - r/c) 



(3.26) 



There is no contribution coming from left side of 1)3.211) thanks to the result 

proved in Eq. ()2.19)1 in the case where we are using a mass-centred frame. The result ()3.25a)l 
generalizes to 2.5PN order the result (4.7a) of paper I. In a future work we shall investigate 
more systematically the relations linking the exterior moments and Sl to the source 
moments II and Jl as they have been defined here. 



IV. EXPLICIT EXPRESSIONS OF THE MULTIPOLE MOMENTS 

Two things must be done in order to obtain the expressions of the moments and Sl 
of Eqs. ()3.25|1 . Firstly one must expand when c — > +oo the source moments II and Jl of 
Eqs. (I3.22|) up to consistent order. Secondly one must evaluate the 2.5PN modification 5Il 
entering the mass moments ()3.25ap . 



A. The source multipole moments 



The post-Newtonian expansion of the source moments 



formed using a formula which was given in Eq. (B.14) of Ref. j29|, namely 



is straighforwardly per- 



^dz5i{z)i:{^,t + z\x\/c) = S(x,t) + 2c2(^2i + 3) 



+ 



5'S(x,t) 



9fS(x,t) + 0(6) . (4.1) 



8c4(2^ + 3)(2£ + 5) 

Explicit expressions of the densities S, Sj and Sjj have also to be inserted into II and Jl] 
these are easily evaluated with the help of Eqs. ()2.8|) . ()2.9|) . ()3.19|) and ()3.23|1 with the result 
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(identical to paper I) 
S = 



C2 



d,Vd,V 



rr 



A/^V - 2VAd,V - W„dlV - -(d,Vf 



+2diVjdjVi + 2d,VdiWj, - -VdiVdiV } , 



(4.2a) 



Si 



' AV 
1 + ^ 



+ 



dkV{diVu - di^Vi) + -dtVd.V \ + 0(4) , (4.2b) 



(4.2c) 



d^Vd.V - -dijdkVdkV \ + 0(2) . 



Note the important fact that the remainder in Eq. 1)4.11) is 0(6) and not 0(5), and thus 
does not contribute to the 2.5PN order (similarly the other remainders in ()4.2|) will not 
contribute). With Eqs. dHH) and (jOJ we recover the same expressions as in Eqs. (4.12) and 
(4.13) of paper I for the source moments in raw form : 



h{t) = FPb=o / rf'x|x 
1 



B 



1 



ttGc^ 



ttGc^ 



XL 



-VdfV - 2VAdiV - W^.dfjV - -idtV) 
+2diV,djV, + 2d,VdiW,j - ^Vd^VdiV 



|X| Xl q2 



2c2(2£ + 3) 



1 4\/\ 
1 + ^ 



I |4 " 



8c4(2£ + 3)(2£ + 5) * 
4{2e+l)xiL 



dfa- 



2(2£+ l)|x|2xi 



c^{£ + l){2i + 3){2i + 5) 



dfa. 



t ^« 



c2(£+l)(2£ + 3) 
2{2i+l)xijL 



c^{e+l){i + 2){2e + 5) 



1 + 



d 



AV 



TVGC^ 



dkV{d,Vk - dkVi) + -dtVdiV 



(7 A A H 

AttG 



diVdjV 



+ 0(6) 



(4.3a) 



Ji(t) = FPB=oeab<it / rf3x|x|-^i XL-l>a ( 1 + ) (Tfo 



2c2(2£ + 3) ^'^^ 



{2i+l)xL-l>ac 

'c^{i + 2){2e + 3) 



dkV{d,Vk - dkVf,) + -dtVdkV 



O'bc + 



AjiG 



-duVdrV 



+ 0(4) . 



(4.3b) 



The remainders 0(6) and 0(4) are negligible. The retarded potentials V , Vi, Wij are then 
replaced by their post-Newtonian expansions when c —>■ +00. It is easily seen that the 
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accuracy of the expansions of V and Wij given in paper I is not sufficient and has to be 
pushed one order farther. The relevant expansions are 



Vi = U, + 0(2) , 



+ 0(2) , 



(4.4a) 
(4.4b) 
(4.4c) 



where the Newtonian-hke potentials f/, X, f/j and P,,- are defined as in paper I by 



[/(x,t) = G / -—-a{^,t) , 
J \^ X I 

X(x,t) = G j rfV|x-xV(x',t) 

f/,(x,t) = G / -^—^a,{^',t) 
/■ rf-'^x' 

P.,(x,t) = G 



X — X' 



cTij + -j—^ (diUdjU - \6ijdkUdkU 



(4.5a) 
(4.5b) 
(4.5c) 

(x,t), (4.5d) 



and where the new terms involve the tracefree quadrupole moment Qij and moment of 
inertia K associated with the mass distribution cr, namely 



(i^xcr(x, t) X. 



K{t) = j d^xa{x,t)x^ . 



(4.6a) 
(4.6b) 



When substituting the expansions ()4.4|) into the source moments ()4.3p all the terms coming 
from the Newtonian-like potentials U, X, Ui and Pij lead to the same expressions as in 
paper I, while the terms coming from the moments Qij and K lead to some correction terms 
(in II only). Let us write 



h = Il + SIl + 0(6) , 
Jl = Jl + 0{4) , 



(4.7a) 
(4.7b) 



where II and Jl are the 2PN-accurate moments which were obtained in paper I, and where 
6Il denotes a 2.5PN correction term (which is distinct from 6Il found in Eq. ()3.25a|) ). After 
the transformation of II as in section IV.B of paper I one can write II and Jl in the form 

4 



hit) 



FP 



B=0 



d xlx 



B 



XL 



cr 



IxPxl 



2c2(2£ + 3) * 
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4(2£ + l)x,x ^ 
"c2(£+l)(2£ + 3) * 

I X I 2^2^ rj4 

'8c4(2^ + 3)(2£ + 5) 

2(2£+ 1) 
'c4(£+ l)(£ + 2)(2^ + 5) 



. 4?7\ 1 
1 + ^ + 



nGc^ \ 4 



2(2£+l)|xpx,^ 
c4(^+ l)(2£ + 3)(2£ + 5) 
1 



t 



(7 a + 

AttG 



diUdjU 



ttGc^ 



-P;,dlU - 2U,mU + 2d,U,d,U, - -{dtUf - Ud^U 



(4.8a) 



JL{t) = FPB=Oeab<i, / d^^\^\^\ XL-l>a ll + ^u]crb + 



2c2(2£ + 3) 



+ 



1 



{2i + 1)XL-I>ac 

"c2(£ + 2)(2£ + 3) 



dkU{dhUk - dM + -dtUdbU 



dt 



CTbc + 



AtiG 



dbUd.U 



(4.8b) 



The moments II and Jl constituted the central resuh of paper I and they were at the basis 
of the apphcation to inspiralhng compact binaries in paper II. On the other hand the 2.5PN 
correction term 61 l is obtained by simple inspection of Eq. (j4.3ap . A simplifying fact in 
obtaining SIl is that the moments Qij and K are only functions of time so that their spatial 
gradients vanish. We obtain 

STl = FPb.o / d'xIxl^l^A-'^rf, + J^Q-f i,a|f/} . (4.9) 

The second term is an integral having a priori a non-compact support. However it can be 
transformed into a manifestly compact-support form by means of the formula 



FP 



B=0 



d^xlxl^- 



2tt 



|x- y| 



2£ + 3 



(4.10) 



This formula is proved by noticing that the integral defined hj Ib = J (i^x|x|'^Xi^|x 



1-1 



is proportional to a prefactor |y| and satisfies AyJ^ = — 47r|y| y^. These two facts imply 



that Ib = Ay ^(-47r|y|^yL) where is defined as in Eq. (3.9) of Ref. j27|, and thus 
that Ib = -47r|y|^+2yi/(5 + 2){B + 2£ + 3) which yields (ICTl after taking the finite part. 
The more complicated formula (4.23) in paper I can interestingly be compared with ()4.10p . 
Thanks to ()4.1()|1 we can write 



Expanding the spatial derivatives in the second term yields finally 



2i + 3 



QfliiQL~i>k - 



2i 



-Q 



(3) 



Kr 



L-2> 



(4.11) 



(4.12) 
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where we have posed 

Qi(t) = j d^xaxL , (4.13a) 
KL-2{t) = J . (4.13b) 

These definitions are in conformity with the earher definitions ()4.6|) . [The brackets <> 
denote the STF projection.] 



B. The 2.5PN modification of the mass moments 

In addition to the previous contribution 51 ^ which is part of the source multipole moments 
Ji, we have seen in Eq. fl3.25a|) that there exists also a 2.5PN contribution 5Il entering Ml. 
Evidently the contribution 5Il is as important as 5Il in that it contributes also to the 
asymptotic waveform depending on the moments Ml. [Note that the terminology refering 
to II and Jl as the multipole moments of the source by contrast to Ml and Sl which are 
viewed as some intermediate moments devoided of direct physical signification is somewhat 
arbitrary. What only matters in fine is to express (by any convenient means) the asymptotic 
waveform in terms of physical quantities belonging to the source.] 

The 2.5PN term 51 l results from the multipole decomposition of the quantity X^^ which 
is itself determined from the other quantity See Eqs. fj3.26p . ()3.14a|) and ()3.11ap . The 
first step in the computation of 5Il is to find the vector ^p^ of the coordinate transformation 
between the inner and external metrics. Since ip^ is of order 0(3,4) and is d'Alembertian- 
free to order 0(7,8) [see Eqs. ()3.2p and ()3.3|) ]. there exist four sets of STF tensors M/l(^), 
XLit), YLit) and ^L(t) such that 



4G 



i>0 



-d, 



-a 



iL 



-WL{t - r/c) 



-XLit-r/c) 
r 



+ 0(7) 



(4.14a) 



+ - 



4G 



E 

i>i 



+ 1 



Ol-i 

£iabdaL-l 



-YiL-i{t-r/c) 



-i(t-r/c) 



0(8) 



(4.14b) 



The powers of 1/c in front of these terms are such that Wl-, ■ ■ ■ iZl have a non-zero limit 
when c +oo. We compute these tensors to the lowest order in 1/c. To do this let us recall 
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the relations between V^*^^* and W^j^*' and the multipole expansions A^(l^) and Ai(Wij) as 
obtained in Eqs. ()3.6b|) and ()3.6cp . We have 



ext 



v: 



ext 



£>0 
£>0 



-dr 



-Vtit-r/c) 
r 



(4.15a) 



-dr 



where we have transformed the relation for W^j into a simpler relation for Vfj, and where 
the moments take to lowest order the form fsee Eqs. (jObjl and ^iE^] : 



rf3xxL(T,(x,t) + 0(2) , 



vf(t) 



Wm = FFb=o / rf^x|xrxiS,,(x,t) + 0(2) . 



(4.16a) 
(4.16b) 



The (non-compact-supported) stress density Sj^ is defined by Eq. fl4.2cp . Having written 
fl4.15|l - ()4.16|l and knowing the explicit multipole decompositions of Vf^^ and V^^^ given 
by Eqs. ()2.11b|l and ()2.11c|) . it is a simple matter to compute the tensors Wl,---,Zl by 
decomposition of the integrands x^cTj and xi'^ij entering ()4.16p into irreducible tensorial 
pieces with respect to their i + 1 and i + 2 indices. We do not detail this computation but 
simply give the result, which is 

2e + i 



Yl 
Zl 



(£+l)(2£ + 3) 
2£ + 1 



d^^XiLai + 0{2) 



2(£+l)(£ + 2)(2£ + 5) 

(£+l)(2£-F3) 
2^ + 1 



FPb=o / rf^x|x|^Xi,xS.,- + 0(2) 

1 



FPB=n / d^y.\y.\^£ab<iiXL~i>bJlac + 0{2) . 



(4.17a) 
(4.17b) 
0(2), (4.17c) 
(4.17d) 



(£ + 2)(2£ + 3)" 

The tensor Wl is manifestly of compact-supported form. There is agreement for this tensor 
with the previous result obtained in Eqs. (2.22a) and (2.19c) of Ref. 29[. As they are 
written the other tensors Xl, Yl-, Zl do not have a compact-supported form. However Yl 
can be re-written equivalently in such a form: 



Yr, 



3 J 3(2f+l) . / 1 

d X I (^^i)(2£+3) ""^*"- - - 



0(2) 



(4.17e) 



The transformation of Yl into the form ()4.17e|) is done using the results (4.2) and (4.18) of 
paper I. 
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The quantities $1"° and X''" can now be evaluated. As is clear from its structure and the 
form of ip^ and M.{y), M-iVi), the quantity VL^^ is made up of a sum of quadratic products 
of retarded waves. We shall write 



00 



-F{t-r/c) 
r 



dr 



-G{t - r/c) 
r 



+ 0(10) 



(4.18) 



where F and G are some functions of the retarded time symbolizing some derivatives of 
the functions ^ V[ and Wl, Xl, Yl, Zl (all indices suppressed). We assume (as can 
always be done) that the derivative operators are trace-free: dp = d^i^di^ ■ ■ -di^^ and dq = 
d^j-^dj^ ■ ■ ■ dj^y. The power of 1/c in front indicates the true order of magnitude of when 
c +00 [see (3.12)], and the remainder 0(10) comes from the un-controlled remainder terms 
in ()4.14|) . To evaluate we need to know the action of the operator FFB=oO]^^r^ni — 1, 
where 1 denotes the unit operator, on the generic term composing Q^^. Actually we shall 
be interested only in that part of X^^ which is strictly larger than the remainder 0(8) we 
neglect in (|3.16|) . In that case, a useful formula shows that all the terms in fi*"^ which are 
composed of the product of two waves with multipolarities p > 1 and q > I yield negligible 
terms in This formula, which is proved in Appendix A, reads 



(FPB=onRV^n - 1) 



Of 











-F(t - r/c) 




-G(t - r/c) 


] 


r 




r 





PQ 



0(3) 



(4.19) 



{Sp^q is the Kronecker symbol). When both p > 1 and q > 1 the right side of ()4.19|) is 
of order 0(3) relatively to the left side and the corresponding term in X^^ is of negligible 
order 0(9). We can thus limit our consideration to the terms in involving at least one 
monopolar wave p = or g = 0. We insert into Eq. ()3.11ap the multipole expansions ()3.7a|) 
and ()3.7b|) together with those of ip^,ip^ given by ()4.14|) . By straightforward application of 
()4.19|) to each of the resulting terms one finds 



X 



00 



16G^>^ 



-dr 



0(9) 



(4.20) 



where the functions W{t) and Yi{t) are given by (j4.17ap in which i = and by (j4.17cp in 
which i = 1, and where we have used the law of conservation of mass implying V^^^ = 0(2) 
and our assumption of mass-centered frame implying Vj-^"* = 0(2). [See the definitions of 
the functions and in Eqs. pi?ajl and ()3.8b;i ; in we denote Vl = and 
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Vl-1 = ^ for the function ()3.8a|) although this notation is shghtly ambiguous with 
fl3.8b|) .] To lowest order the function Viit) reduces to 

Vl = Ql + 0{2), (4.21) 

where Qiit) is the moment defined in Eq. ()4.13a|) . On the other hand W{t) satisfies 

W=^J d'^Xia, + 0(2) = i + 0(2) , (4.22) 

where K[t) is the moment of inertia ()4.6b|) . Similarly one finds 

Y, = ^Gf^ +0{2) (4.23a) 

where the vector Giit) reads 

Oj = y (fyi ^cTiX^ - ^a-jXiXj^ . (4.23b) 
With these notations we end up with the 2.5PN correction term 5Il [compare ()3.26p and 

(11201)], 

5h = § {Ik^'^Ql - iKi'^Qi' + f <Ql-i>} + 0(7) . (4.24) 

Summarizing the results so far, we have explicitly computed the mass multipole moment 
Ml given by Eq. ()3.25a|l . It contains a 2.5PN contribution issued from the source moment II 
and computed in Eq. ()4.12|1 . and also the direct 2.5PN modification computed in Eq. ()4.24j) . 
We can write 

Ml = Il + A/l + 0(6) (4.25) 

where II is given by Eq. ()4.8aj) (this was the result of paper I), and where AIl = SIl + SIl 
is given by 

-^Oa.«t-i>* + |Gg,«L-i>} + 0(7) . (4.26) 

We recall that the tensors Ql, Kl_2 and Gi are defined in Eqs. Kl^ and (103^)) . The low 
orders in i read 

AI = —MK^^^ + 0(7) , (4.27a) 
3c^ 

A/, = l^MOf +0(7) , (4.27b) 
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A/.- = § [^K^'^Q^, - ^K^Q^^ - ^KQ^ - ^Qag,>.} + 0(7) , (4.27c) 

where M is the total mass such that Q = M + 0{2), and where we have used a frame such 
that Qi = 0(2). The quadrupolar correction term ()4.27c|) will contribute to the asymptotic 
waveform at the 2.5PN order. The dipolar correction term (j4.27b|) will be used to determine 
the center of mass of the system at this order. 



V. THE 2.5PN-ACCURATE GRAVITATIONAL LUMINOSITY 



The mass and current multipole moments Ml and Sl are determined up to the neglect of 
0(6) and 0(4) terms respectively and can be used to compute the gravitational luminosity 
(or energy loss rate) of the system at 2.5PN order. To compute the waveform at the same 
order would necessitate a more accurate determination of the current moments, up to the 
neglect of 0(5) terms. We shall leave this computation for future work. 

Let = {cT, X) be a coordinate system valid in a neighbourhood of future null infinity 
and such that the metric admits a Bondi-type expansion when R = |X| +oo with 



Tr = T — R/c staying constant. See j42| for the proof (within the present formalism) of the 
existence and construction of such a coordinate system. The relation between Tr and the 
retarded time of the harmonic coordinates xfi^ is 



can 

Tn = tcan - ^ - ^ In (^) + 0(l/r_) + 0(5) , (5.1) 

where b is some arbitrary constant time-scale. [Actually, to be consistent with the 2.5PN 
precision one should consider also the next-order post-Newtonian term in Eq. ()5.ip ; however 
this term will not be needed in the following.] It is sufficient to control the transverse and 
tracefree projection of the leading-order term ~ R~^ in the spatial metric. A multipole 
decomposition yields a parametrization into two and only two sets of STF moments Ul and 
Vl which depend on and can be referred to as the "radiative" or "observable" mass and 
current moments. These are chosen so that they reduce in the limit c +oo to the ith 
time derivatives of the ordinary Newtonian mass and current moments of the source. The 
total luminosity C = C{Tji) of the gravitational wave emission when expressed in terms of 
Ul and Vl reads 



+ 00 



1=2 



c2^+i \(£- 1)££!(2£+ 1)!!^ ^ {i - + l)\{2i + ly.lc^^ ^ I • 
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Considering L to 2.5PN order one retains in ()5.2j) all the terms up to the neglect of a 
remainder 0(6), and finds 



1 

+ ^ 



' rr(l) 7-/(1) I ^yy'iyy 
9072 34 u'fc ^i'^ 



1 



{1)t/(1) 



+ 0(6) 



(5.3) 



Because the powers of 1/c go by steps of two in £, this expression is in fact the same as 
already used in paper I. 

All the problem is to find the relations between the radiative moments Vj, and the 



moments M^,, 5*^ we have previously determined. We rely on previous papers 
paper I) having written the general form of these relations as 
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G 



n— 1 



n>2 



I 1 G^~^ 



n>2 



g3(n-l)+S£i-£ 



and 

(5.4a) 
(5.4b) 



where the functions XnL and Yni represent some nonlinear (and in general non-local) func- 
tional of the moments Ml and Sl- The powers of 1/c in Eqs. fj5.4|) come from the di- 
mensionality of the functionals XnL and YnL which is chosen to be that of a product of n 
multipole moments and their time derivatives. We can write symbolically 



{a„) 



(5.5) 



The notation in Eqs. (j5.4j) - (j5.5p is the same as in paper I ; in particular S£j denotes the 
total number of indices on the n moments in the term in question, and denotes the sum 
S£j + s where s is the number of current moments. Here we shall need only the fact that 
S£j is larger than the multipolarity i by an even positive integer 2k which represents the 
number of contracted indices between the moments composing the term (with the current 
moments carrying their associated Levi-Civita symbol), i.e. 



+ 2k . 



(5.6) 



With the latter equation it is simple to control the type of nonlinearities which are present in 
the radiative moments to 2.5PN order. Since the reasoning has already been done in paper I 
to 2PN order we consider only the case which is further needed, that of a term of pure 2.5PN 
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order in the mass- type quadrupole moment Uij (having £ = 2). By Eqs. ()5.4ap and ()5.6p this 
case corresponds to 3{n — 1) + 2A; = 5. The only solution is = 2 (quadratic nonlinearity) 
and k = 1 (one contraction of indices between the moments). With two moments, one 
contraction and i = 2 one has ii + i2 = 4. Furthermore one of the two moments is non- 
static, hence 1,2 > 2 say, so we obtain only two possibilities {i-^^, £2) = (1, 3) or (2,2). The first 
possibility is excluded because the moment having i^ = lis necessarily the constant mass 
dipole Mi which has been set to zero. Only remains the second possibility li = i.2 — which 
corresponds either to the interaction between two mass-type quadrupole moments Mij or to 
the interaction of the (constant) current-type dipole Si with Mij. 

We combine these facts with the study done in Ref. 3^ of the occurence of "hereditary" 
terms in the asymptotic metric at the quadratic approximation n = 2. Two and only two 
types of hereditary terms were found: the "tail" terms coming from the interaction between 
the monopole M and non-static multipoles, and the nonlinear "memory" term which is made 
of the interaction between two non-static multipoles. The tail terms are of order and 
have been included in paper I, but the memory term arises at the order (in the radiative 
quadrupole Uij). The latter term can be straightforwardly computed from Eqs. (2.42a), 
(2.21) and (2.11a) in Ref . H- 

An equivalent result can be found in Ref. [4^. For discussions 
on the memory term see Refs. It is clear by the previous reasoning that the memory 

term represents the hereditary part of the 2.5PN contribution in the radiative moment Uij, 
corresponding to the interaction of two moments Mij. Associated with this term there 
are also some instantaneous terms having the same structure and exhausting (a priori) the 
possibilities of sharing time- derivatives between the two moments. 



Gathering these results with the results of Ref. 



and paper I we obtain the expression 



of the radiative quadrupole Uij to 2.5PN order as follows : 
f/,(T,) = <)(T,) + ^ £^ dr [in Q + - 

+/3Mi%M^'J, + iMflM,^, + A^,Mit,e,>,^| + 0(6) . (5.7) 

The constant h in the term of order (tail integral) is the same as in Eq. 1)5.111 . The memory 
term is the integral in the brackets of order c~^. The qualitatively different nature of these 
two integrals can be clearly understood when taking the limit Tr +00 corresponding to 
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very late times after the system has ceased to emit radiation. In this hmit the third and 
higher time derivatives of Mij are expected to tend to zero, so the tail term tends to zero 
while by contrast the memory term tends to the finite limit -{2G/7c^) J^^ duMl%M^% 
(see Refs. 



45l|). The coefficients a, (3, 7 and A are some purely numerical coefficients in 
front of instantaneous terms (which depend on only). These coefficients can be obtained 
by a long computation using the algorithm of Ref. j27|, however we shall not need them in 
the application below (they will be computed in a future work). The higher-order radiative 
moments take similar expressions but are needed only to a lower precision. The relevant 
expressions for Uijk and Vij have been written in paper I, and read 



2GM 



dr 



dr 





97" 


Ki) 






^ 60_ 



Mlf,{Tn-r) + 0{b) , (5.8a) 



c JO 

while the other needed moments are given by 





7" 








+ 6. 



sI^\Tr-t) + 0{5) , (5.8b) 



UijkmiTji) 

VijkiTii) 



Mi;ijTn) + 0i3) 
S§liTR) + 0{3) . 



(5.9a) 
(5.9b) 



The expressions fl5.7j) - ()5.9|l of the radiative moments are to be inserted into the gravi- 
tational luminosity ()5.3|) . This leads to a natural (though not unique) decomposition of C 
into instantaneous and tail contributions, 



inst ~r ■'-tail 



(5.10) 



The contribution Ci^ist depends only on the instant Tr and is given by 



r, _ G (I (3) ^(3) , 1 r 1 yM-(4) y^^(4) 16 ^(3) ^(3) 

~ '^{5 189 '^'^ 45 

5c5 



H M^'- 



[a 



^/M-(5) »^(5) 1 ^(4)^(4) 

9072 ^^^^ ^^^^ 84 ^"'^ 



HP + i)Ml^Mf^ + ^M^l^M.k +XSkM^^,e,kn\ + 0(6)| . (5.11) 

Note that this involves the term coming from the nonlinear memory which is instantaneous 
in the energy loss (the memory effect exists only in the waveform). The tail contribution 
depends on all instants Tr — r anterior to Tr and reads 
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+i^S.?'(r«)/"°°<i-S.?'(rH--)ln(^) +0(6)} . (5.12) 
where we have posed for simphcity 

6l = 6e-"/l^ 62 = &e-9^/6°, 63 = 6e-^/^ (5.13) 

The luminosity ()5.1()|1 - ()5.13|1 in which the moments have been determined in Sect. IV is our 
final result for the general case of a (semi relatistic) isolated system. 



VI. APPLICATION TO INSPIRALLING COMPACT BINARIES 

The authors of Ref. j2l| (paper II) applied the results of paper I to an inspiralling compact 
binary system modelled by two point-masses moving on a circular orbit. Here we do the 
same for the results derived previously and obtain the energy loss rate and associated laws 
of variation of the frequency and phase of the binary to 2.5PN order. 



A. The equations of motion 

Obviously the equations of motion of two point-masses at the 2.5PN approximation are 
needed for this application. These equations have been obtained in the same coordinates 

n 

as used here by Damour and Deruelle [S2| studying the dynamics of the binary pulsar. For 
inspiralling compact binaries one needs only to specialize these equations to the case of an 
orbit which is circular (apart from the gradual inspiral). A relevant summary of the Damour- 
Deruelle equations of motion is presented for the reader's convenience in Appendix B. Here 
we quote the results valid for circular orbits, following mostly the notation of paper II. Mass 
parameters are denoted by 

m = mi -f m2 ; Xi = — ; X2 = — ; v = X1X2 (6.1) 

m m 

(the total mass is henceforth denoted by m to conform with paper II). The individual 
positions of the two bodies in harmonic coordinates are y\ and 7/2- Their relative separation 
and relative velocity are 

x' = y\-y\ ; = ^ ■ (6-2) 
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A small ordering post-Newtonian parameter is defined to be 

7 = ^ 6.3 

with r = |x|. Next we use the fact that the origin of the coordinate system is located at the 
center of mass of the binary. This means that Mj = where Mj is the dipole mass moment of 
the external field. By Eqs. ()4.25|1 and ()4.27bj) this means Jj+AJj = 0(6) where Jj is the dipole 
moment which was computed in paper II (before expressing it in the relative frame) and 
where AJj = (4G/5c^)mG'^^'' + 0(7) where Gi = J d'^x (o"jX^ — ^ajXiXj). For circular orbits 
one finds = mi/(X2-Xi)rV+0(2) and thus A/^ = (AG /5c^)Gm^u{Xi-X2)ivyr)+0{7). 
This readily shows how the relations (3.7) of paper II are to be extended to 2.5PN order. 
We find 

yl = [X2 + Si^^'iX, - X2)]x' - -^^^^(Xi - X2y + 0(6) , (6.4a) 

5 rc^ 

4 G'^m^u 

= [-X, + 3z/f (Xi - X2)]x' - -^^^{X, - X^y + 0(6) . (6.4b) 

TC 

This result is in agreement with the 2.5PN-accurate center of mass theorem of Refs. |37|. 
The assumption that the orbit is circular apart from the adiabatic inspiral due to reaction 
effects of order 0(5) implies that the scalar product of x^ and f* is of small 0(5) order: 
x.v = [xv) = 0(5). By Eqs. ()6.4j) this implies also (nvi) = 0(5) and {nv2) = 0(5) where 
ra* = x^/r. These facts simplify drastically the equations of motion of the binary given by 
Eqs. (B1)-(B2) in Appendix B. The result when expressed in the relative frame simply reads 

^2PN^ - — — + 0(6) , (6.5) 



dt 5c^ r 

where we have introduced the angular frequency co'2pn defined by 



2 _ Gm 



1- (3-z/)7 + ( 6 + ^z/ + z/2 ) 7^ 



(6.6) 



This frequency represents the orbital frequency of the exact circular periodic orbit at the 
2PN order (see Eq. (3.11) in paper II). The relation between the norm of the relative velocity 
V = |v| and u;2pn is obtained by multiplying both sides of ()6.5|) by x\ Using (xv) = 0(5) and 
d{xv)/dt = 0(10) (because d{xv)/dt is of the same order as the square of reaction effects) 
we find 

w = ra;2PN + 0(6) . (6.7) 
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Finally we write the result for the orbital energy E = £;2.5PN g^tering the left side of the 
energy balance equation dE/dt = —C^ + 0(6) derived in Eq. ()B10|) of Appendix B. This 
energy is computed from Eqs. ()B6|) and ()B11|) in which one uses the circular orbit assumption 
together with (j6.7j) . The result is 

E = -—mu-f j 1 - - - g(7 - 49i^ - '^'^h^ > + 0(6) . (6.8) 

There are no terms of order 7^/^ = 0(5) for circular orbits because the term 0(5) in Eq. (Bll) 
is proportional to (nv) and thus vanishes in this case. 



B. The energy loss rate 

The gravitational luminosity £ of a general source was split into two contributions, an 
instantaneous one £inst given by Eq. ()5.11|1 and a tail one £taii given by Eq. ()5.12j) . We shall 
basically show that only £taii contributes to the 2.5PN order in the case (but only in this 
case) of a binary system moving on a circular oihit. 

Let us consider first the contribution Ci^st- The only moment it contains which is required 
with full 2.5PN accuracy is the mass quadrupole moment 

Mij = 4- + A/,, + 0(6) (6.9) 

where lij results from Eq. ()4.8ap and A/jj is given by Eq. ()4.27cp . For a circular orbit, the 
moment of inertia K is constant [indeed K = mz/r^ + 0(2)], hence Alij reduces to two terms 
only: 

A/., = § - lQfl^Q,>k] + 0(7) . (6.10) 

We prove that the contribution in £inst which is due to AJjj is in fact zero. Indeed this 
contribution is made of the contracted product between qI^^ and Alij^ (recall that M^- = 
Qij + 0{2)), hence of the contracted products Qf^Qf^ and Qfj^Q% '^^Q^/k- But in the 
circular case an odd number of time- derivatives of Qij yields a term proportional to x^'^v^^ 
while an even number yields either x'^^x^^ or v'^^v^^ . Thus a contracted product of the type 
Qij'^Qif where n + p is odd necessarily involves one scalar product (xv) and is thus zero; 
similarly a product of the type Qij'' Qik' Q^jk where n + p + q is odd is also zero. The same is 
true of a product like ^ijkQ^mQkm where r + s is even. These simple facts show that ()6.10|) 
cannot contribute to the energy loss rate. Furthermore we find that all the terms in £inst 
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which involve the contracted products of three moments [i.e., all the terms of order c ^ in 
Eq. ()5.11|) ] are also zero. Hence we can write in the circular orbit case 



G n 



inst 



- r{4) 7(4) 21. n-'' T 

Igg ijk ijk ' 45 ' 



16 



r(3) 7(3) 



r(5) 7(5) 



'(4) 7(4) 



9072 '■''^"^ ^-''^"^ 34 «i 



+ 0(6) 



(6.11) 



Now recall that the moments 1^ and are the ones which were used as starting point in 
the computation of paper II. In using paper II one must be careful that in this paper the 
equations of motion have a precision limited to 2PN instead of 2.5PN. We first note that 
the 2.5PN terms in the mass-centered-frame relations ()6.4|1 will cancel out when expressing 
the quadrupole mass moment in terms of the relative variables. Furthermore during the 
computation of this moment in paper II the equations of motion were used only to reduce 
some terms which were already of IPN order. Thus all the expressions of the moments 
II and Jl computed in paper II (but not the expressions of their time derivatives) can be 
used in the present paper without modification. Notably the mass quadrupole lij takes the 
expression (3.74) of paper II, modulo negligible 0(6) terms: 

„2 



7 



7 



— l + 39z/ - 
42^ ' 1512 



(461 + 18395Z/ + 241z/' 



11 



7 



2r'-''^^ + 378 



(1607 - 1681z/ + 229i/^ 



+ 0(6) .(6.12) 



However we need the third time derivative of ()6.12j) which does involve, because of the more 
precise equations of motion, some new terms with respect to paper II. We find 



r(3) 



STFjj vm{ — < 



Gm 



1 - ^(149 - 69z/) + ^—(7043 - 7837z/ + 3703z/^ 
42 1512 



+ 



64 G^m^u 



Gm 



0(6) . 



(6.13) 



But for the same reason as before these new terms will not contribute to the energy loss 
for circular orbits. Thus we conclude that £inst is exactly given by the expression found in 
Eq. (4.12) of paper II modulo negligible terms of relative order 0(6) = 0(7^). We can thus 
write 

32c5 



-"inst 



5G 



■z/V<! 1 



/2927 5 
V^36" ^A^^^ 



293383 
9072 



^^7^ + 0(7^^ 



(6.14) 
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Let us now turn our attention to the tail part of the gravitational luminosity as defined 
by Eq. ()5.12p . The correct formula to compute the tail integrals in £taii is 

^'°°drl„(^)cos(nr) = -|^. (6,15) 



where f2 denotes the (real) angular frequency of the radiation. It was proved in Ref. 
that this formula is to be applied as it stands (i.e., even though the integral is not absolutely 
convergent) to a fixed (non-decaying) orbit whose frequency is equal to the current value 
of the orbital frequency uj = ci;2Pn(7j?). The numerical errors done in applying this formula 
were shown to be of relative order 0(lnc/c^), which is always negligible for our purpose here 
(because of the explicit powers of c^^ in front of the tail integrals). Thus we replace into 
the integrand of ()5.12|) the moments Mi and Sl by their expressions valid for fixed circular 
orbits and up to the appropriate precision. The necessary formulas (taken from paper II) 
are 



< =32STF,,^xV 



1 - ^(149 -69z/) +0(4), (6.16a) 



42 
42' 



1 - — (275- lllz/) +0(4), (6.16b) 



M^i = 3(X2 - Xi) STF,,fc {7^^*'' - 20xV'=| + 0(2) , (6.16c) 

3(X2 - Xi) STFi.-fc — ^ <^ -Qi—-x'jk + l82xV'= \ + 0(2) , (6.16d) 



= -(X2 - Xi) STFi,- -^e-^^v^x'^v'' + 0(2) , (6.16e) 



Grn^v 

■i2 3 

Sf = (X2 - Xi) STF,, ^^-^£'^%^x''v' + 0(2) . (6.16f) 



The quadrupole moment includes IPN terms while the other moments are Newtonian. The 
computation of £taii involves many scalar products between or f * given at the current 
time Tr, i.e. = x*(Tr) or = ^^(Tr), with the same vectors but evaluated at an arbitrary 
earlier time Tr — r, say x'* = x*(Tr — r) and f'* = f*(TR — r). Relevant formulas for these 
scalar products are 

{xx') = r^cos(u;r) , (6.17a) 
{vv') = r^cj^ cos(u;r) , (6.17b) 
{xv') = —{vx')=r'^ujsm{ujT) , (6.17c) 
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where u is the current value of the frequency. The reduction of £taii is quite straightforward. 
We need to use uj"^ = {G'm/r^)[l — (3 — z/)7 + 0(7^)] (the IPN correction in u;^ is sufficient) 
and the formula (|6.15|) where fl = u, 2lj or 3u!. One uses also {X2 — Xi)^ = 1 — iu. The 
result reads 

32c^ 2 , r 0/2 /25663 125 \ ,,2 ^/ sxl . ^ 

Aaii = ^^ 7' |47r7^/' - + — 7r7^/' + 0(7=^) | . (6.18) 

The complete 2.5PN-accurate gravitational luminosity generated by a inspiralling com- 
pact binary moving on a quasi-circular orbit is therefore obtained as the sum of Eqs. ()6.14|1 
and (IFTTHIl . We obtain 



32c^ 2 J /^2927 5 \ 0/2 /293383 380 \ 

5G ' \ V 336 4 y ^ ' V 9072 9 J 



V I 7^ 



/25663 125 ^ 5/2 ^, 3, 1 

+ — z/ 7r7^/^ + 0(7^)L (6.19) 



V 672 

where we recall that the post-Newtonian ordering parameter is 7 = Gmjr(? with r being the 
harmonic-coordinate orbital separation. This expression was already obtained to IPN order 



in Refs. m 
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)ordmate orbital separation, inis expression was already obtamea to iriN order 
42, 3. to 1.5PN order in Refs. 00,0,0. and to 2PN order in Refs. [3,0 



22]. The 2.5PN order added in this paper is like the 1.5PN order due to the presence of 
the radiation tails in the wave zone, as indicated by the irrational number vr in factor coming 
from the formula (6.15). Hence there is up to 2.5PN order a clean separation between the 
integer post-Newtonian approximations which come from instantaneous relativistic effects 
in the source multipole moments and the half-integer approximations which are due to 
hereditary effects in the wave zone. However the next post-Newtonian approximation (3PN) 
is expected to involve both types of effects. 



C. The orbital phase 



Once the energy loss has been derived in a particular coordinate system it can be re- 
expressed in a coordinate-independent way by using the directly observable frequency lo = 
Ci;2PN instead of the parameter 7. Defining x = (Gmcu/c^)^/^, we obtain the inverse of 
Eq. dnini) as 



7 



X 



1+ I 1-- lx + 



65z/ 



which is substituted into Eq. (j6.19j) with the result 

32c^ 2 J /1247 35 \ o,, 
^'^^^- ^+12'^'^ + ^"^' 



+ O(x^) 



(6.20) 



44711 9271 65 2 \ 2 
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An important check is obtained in the test-body hmit z/ ^ of this resuh which is found 
to agree with the 2.5PN truncation of the resuh of perturbation theory known through the 



equivalent of 4PN order 15|, [l^ llZl (see the equation (43) of Ref. |17|). As for the orbital 
energy E = E'^'^™ which has been given in Eq. ()6.8p following the equations of motion 
summarized in Appendix B, it reads in terms of X clS 



E 



—mux \ 1 - —(9 + z/)x - - ( 27 - 19z/ + — ) + 0{x 



(6.22) 



The laws of variation of the frequency u and phase of the binary during the orbital 
decay will be computed using the energy balance equation 

dE 



dT, 



R 



(6.23) 



where Tpt is the retarded time in the wave zone. Note that we meet here the remarks made 
in the introduction. Up to now the analysis has been rigorous within the framework of post- 
Newtonian theory. Namely, C has been computed from a well-defined formalism based on 
convergent integrals for the multipole moments (see Section IV) . The reduction of C to binary 
systems using delta functions to describe the compact objects can probably be justified at 



the 2.5PN order by using the results of Ref. 



Ill ] (see paper II). Furthermore E is computed 



directly from the D amour- Deruelle equations of motion • We now postulate the validity of 
the energy balance equation (|6.23|) where both E and C take their full 2.5PN accuracy. This 
equation has been proved here (in Eq. ()B9|) of Appendix B) only when C takes its ^ 



value C . It was proved before to Newtonian order in Refs. 



33, 
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other hand the balance equation is also known to hold at IPN order 



tails at 1.5PN order m 
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36, 



ewtonian 
37|. On the 
and even for 



40|. To prove it at 2.5PN order like in ()6.23p would mean knowing 
all the radiation reaction effects in the binary's equations of motion m 
To palliate this it has been customary in this field (see Refs. 



:o the 5PN order. 
1^ 21|) to 



compute the binary's orbital decay assuming that the energy balance equation is valid to 
high order (the angular momentum balance equation is unnecessary for circular orbits). 

Adopting the same approach to the problem, we introduce like in paper II the adimen- 
sional time 



e 



bGm 
36 



R 



(6.24) 



Equation ()6.23|) is then transformed into the ordinary differential equation 

/34103 13661 59 2' 

+ \ ly H ly I 

V 18144 2016 18 



On the other hand the differential equation for the instantaneous phase (p of the binary is 
d(p = LvdTR = (h/v)x^^'^d6 or equivalently 



dx 64u V336 4 



/3058673 5429 617 

+ \ z/H z/" I 

V 1016064 1008 144 



+ (6.26) 



The frequency and phase follow from the integration of these equations. The frequency (or 
x-parameter) reads 



X 



4 \ V4032 48 y 5 



-3/8 



/ 19583 24401 31 , \ ^ 1/2 

+ \ ly H ly <d ' 

V254016 193536 288 



+ 



/_11891 A , , I 

V 53760 1920 J / 



where = O^. — 9 denotes the (adimensional) time left till the final coalescence, and the 
phase is 



32z/ i V1008 12 

/ 15293365 27145 3085 o \ 1/2 



1016064 1008 144 



/38645 65 \ , /x\ ^. ^i^A ,^ 

where the constant xq is determined by initial conditions. By substituting Eq. ()6.27|1 into 
()6.28|) one obtains the phase as a function of time: 

^ z/\ V8064 96 y 4 
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/ 9275495 284875 1855 2\^i/s 
V 14450688 258048 2048 ) 

( 38645 65 \ , / 6 \ .i^A , , 

where Go denotes some constant. 

The presence of the logarithm of the frequency or coalescing time in the phase is a novel 
feature of the 2.5PN approximation. As a result the phase tends formally to infinity at 
the time of coalescence. In other words, the number of cycles left till the time of coalescence 
starting from an initial frequency is always infinite. This is in contrast with the previous 
2PN approximation where was tending to a finite value (^c- Of course this peculiarity is 
not relevant to the actual physics of the binary because the post-Newtonian approximation 
breaks down before the coalescence, at some finite value of the frequency. The fact that 
tends to infinity is merely a consequence of our idealized model assuming that the two bodies 
have no internal structure or horizons, which is valid only during the inspiral phase preceding 
the coalescence. However, during this phase the model is expected to be very precise and 
in particular the logarithmic 2.5PN term in the phase has to be taken numerically into 
account. This logarithm makes in fact the contribution of the 2.5PN term to the total 
number of cycles in the detector's bandwidth somewhat more important than what could 
be deduced from a simple order of magnitude estimate. 

Let us be more precise and evaluate the contribution of each post-Newtonian terms in 
the phase ()6.28|1 to the accumulated number M of gravitational- wave cycles between some 
initial and final frequencies tUj and lOj. Note that such a computation can only be indicative 
of the relative orders of magnitude of the different terms in the phase. A full analysis 
would require the knowledge of the power spectral density of the noise in a detector, and 
a complete simulation of the parameter estimation using matched filtering P, IS, 1^ 1^ ■ 
The contribution to M due to the 2.5PN term in (16.2811 is 



48z/ V 1344 16 y V^i/ 
In the case of the inspiral of two neutron stars of mass 1.4M0, and with the values cjj/vr = 
lOHz (set by seismic noise) and u;//7r = lOOOHz (set by photon shot noise), we find from 
Ref. |jj and Eq. (On|l 
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Newtonian 



16,050 



IPN 



439 



1.5PN 



-20^ 



2PN 



9 



2.5PN 



-11 



[Note that as far as the final frequency cu/ is determined only by detector's characteristics 
(and not by the location of the innermost stable orbit), the term A/'2.5pn depends on the 
masses only through the mass ratio v.] 

In this (indicative) example we see that the contribution of the 2.5PN order more than 
compensates the contribution of the previous 2PN order. [But of course both contributions 
have to be included in the filters since they have different functional dependences on the 
frequency.] However note that finite mass effects (proportional to v) are numericall y sm all 



at the 2.5PN order, contrarily to the 2PN order where they are quite significant |12,|21|. In 



the previous example, they contribute only to —0.1 cycle as compared to A/'2.5pn = —11- 



APPENDIX A: THE PROOF OF SOME FORMULA 



To prove the formula (j4.19|) we begin with the identity 

r^Uf = □(r^/) - B{B + l)r^-^f - 2Br^-^drf 
which permits the transformation of the left side of ()4.19|) . namely 



pXq = (FPB=onRV^n ~i)(dp 



-F{t-r/c) 
r 



dr 



-G{t - r/c) 
r 



into 



pXq = FPB=on]j' -Br^'^[2dr + r'^] dp 



'-F{t-r/c) 



dg -G{t-r/c) 



(Al) 



(A2) 



(A3) 



Because of the presence of the explicit factor B in the integrand, the equation ()A3|) appears 
to be the residue of the Laurent expansion of the retarded integral near 5 = (we have 
discarded a term proportional to B^ which is known to give zero contribution). Residues 
of retarded integrals have been investigated in Eq. (4.26) of Ref. [4^ and yield retarded 
solutions of the wave equation. By a simple dimensional analysis we find that pXq admits 
the structure 



[(P+'?)/2] 

E 

fe=0 



-F^''\t-rlc)G^^\t-rlc) 
r 



(A4) 
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Here, [{p + q)/2] is the integer part of {p + q)/2, Ol denotes a tracefree derivative operator 
composed of £ spatial derivatives where I = p + q — 2k, and the number of time derivatives 
on F and Gisa + 6 = 2A; + 1. Omitted in (A4) are pure dimensionless coefficients and 
numerous Kronecker symbols. The explicit dependence on 1/c present in Eq. ()A4|) shows 
that if we are interested only in the leading-order term when c — > +00 we can limit us to the 
computation of the terms having the maximum number of space derivatives, i.e. £max = p + q 
(corresponding to k = 0). The errors made in keeping only such terms are of relative order 
0(3). We use to compute ()A3|) the formula 



-Fit-r/c) 
r 



Up 2_ a; 
=0 



F(P-^)(t-r/c) 



(A5) 



where the numerical coefficient is 



(-)P(p + z)! 



(A6) 



2H\{p-i)\ 

(see Eq. (A35a) in [27*]). This yields a double sum of terms having in front the product of 
the tracefree tensors hp and hq. This product can be decomposed into tracefree tensors of 
multipolarities i = p + q — 2k with k = 0, [{p + q)/2], however if we are looking only for 
the term having ^max = P + q one can simply replace the product hphg by the tracefree 
tensor hpq (with one as a coefficient in front). Hence 

F'-P-'^t - r/c)G^''-^\t - r/c) 



pXq = ^afajFPB=on^' <^ -Br''-'hpQ[2dr + 



+ 0(3) , 
(AT) 

where the remainder 0(3) comes from the errors made in the replacement hpfiq ^ hpq. The 
integrand in ()A7|) is easily transformed into a sum of terms of the form r^~^hpQH(t — r/c) 
times the factor B, where K is an integer. Usi ng t he fact that the residues of such terms 
exist only when K > p + q + 3 (see Eq. (4.26) in j40|), one finds that the summation indices 
i and j in ()A7|1 must satisfy i + j = p + q or p + q — 1. This simplifies much the computation 



of 
in 



7j) which is done using the expression (A6) of the coefficients and the equation (4.26) 
4(1. The result is 



d 



PQ 



rc 



- 5,,oF^'^G + 6o.,FG^'^ 



+ 0(3) 



(A8) 



where 6pfi denotes the usual Kronecker symbol. 
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APPENDIX B: THE DAMOUR DERUELLE EQUATIONS OF MOTION 



We present a summary of results concerning the equations of motion of two point-masses 
moving under their mutual gravitational influence up to the 2.5PN order. These equations 
were obtained by Damour and Deruelle js^l (see also the presentation [ll| of which we are 
close here) . The equations of motion of body 1 (say) take the Newtonian-like form 



dt 
dv\ 
'dt 



^i + 3^i + 4c^i + i^i+o(6) 



(Bla) 
(Bib) 



c 

where y\ and v\ denote the instantaneous position and coordinate velocity of body 1 (in the 
harmonic coordinate system), A^^ is the usual Newtonian acceleration of the body, and 
C\ and D\ represent the relativistic corrections of order IPN, 2PN and 2.5PN respectively. 
The equations (12.7)-(12.12) in Ref. [ll| give these terms as 

Gm2 



A\ 
B\ 



Gm-2 



-n 



n 



-v\ - 2vl + ^{vxV2) + ^{nv2f + 5 



Grux Gm 



+ 4- 



+(^i-^2) 



4(nt>i) — 3{nv2) 



(B2a) 



(B2b) 
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Gm2 ( 



1 



n' -2v\ + ^vl{vxV2) - 2(t;if2)^ + -;^v\{nv'^'^ 
^-%v\{nv2f - ^{vxV2){nv2f - ^{nv2f 



Grrii 
r 

Gm2 



17 



+ ^v^ - ^{viV2) + Y(nt;i)^ - 39{nvi){nv2) + y(™2)^ 



4:vl - 8{viV2) + 2{nvif - 4:{nvi){nv2) - 6{nv2y 



vl{nv2) + 4t>2(nt>i) — 5^2(^2) — 4(fif2)(nfi) 



9 



+A{viV2){nv2) - Q{nvi){nv2)'^ + -{nv2Y 



+ 



Grrii 



Grri') 



-2{nvi) — 2{nv2) 



G^m2 



-n 



57 



69 



~—m^ — 9m2 —mim2 



D\ 



4 G'^mim2 
5 



r, Gm^ 
-f 2 + 2 



Gm-i 



+ n\nv) 



' 2 ^Gmi 52Gm2 
3v — 6 h — 



(B2c) 
kp2d) 
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where mi and m2 are the two masses, vl and the two velocities, r = |yi — y2| the 
separation, and = {y{ —y\)/r. Scalar products are denoted by e.g. vr V2 = (^1^2). In the 
last equation t>* = v\ — v\. The equations of the second body are obtained by exchanging 
the labels 1 and 2, being careful to the fact that rt and change sign in the exchange. The 
acceleration term D\ is responsible for the dominant damping or radiation reaction effects 
in the dynamics of the binary. 

The equations of motion when truncated to 2PN order (neglecting the damping term 
D^) admit a Lagrangian formulation and the associated conservations laws [37]. The 2PN 
Lagrangian in harmonic coordinates depends not only on the positions and velocities of the 
two bodies, but also on their accelerations a\ 2 = dv\ 2/dt: L^™ = L^™(?/, t>, a). It is given 
by 



'2PN 



{y,v,a) 



Gm 17712 



1 Gmim2 1 



H — miv-^ 



7 



1 



-fi - -{V1V2) - -{nvi){nv2) - - 



1 Gmi 



+— 7771771 + 

lo 



2 ' 4' ' 4 

Gr77i7772 [7 4 15 2 2 o 2/ N 

7, 



+ o('"i'"2) - -{nvi) V2 + -{nvi){nv2){viV2) + —{nvi) (77772) 



4 



16 



G'^mfm2 



1 7 7 7 1 7 

■^vl + -vl - -{V1V2) + -{nvif + -{nv2f - -{nvi){nv2) 



+Gmim2 
G^mim2 



(nai) ( ^772 - ^(^^^2)^ 



7 

-(772ai) (77772) 



1 2 1 2 19 

-771^ + -7772 + -^rnim2 



(B3) 



The summation symbol runs over the two bodies 1 and 2. The equations obtained by 
variation of this Lagrangian, and in which the accelerations are replaced after variation by 
their (Newtonian) values, are equivalent to the 2PN equations of motion (Eqs. (B1)-(B2) 
where Dl = 0) modulo negligible 0(6) terms. From this Lagrangian one constructs the 2PN 
integral of energy 



E 



2PN 



dL 



2PN 



d dL 



2PN 



+ a\- 



dL 



2PN 



-L 



2PN 



dvl dt da\ J da\ 
and one replaces in it the accelerations by their (Newtonian) values, resulting in 



(B4) 



(B5) 



42 



The explicit expression of the energy E"^^^ as a function of the positions and velocities is 



772PN/ N \^ j 2 lGmim2 3 4 
E {y,v) = 2^^^m,v, - + ^rn,v^ 



+ 



Gmim2 



3 „ 7 
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(B6) 



The fact that E"^™ is the integral of energy of the 2PN equations of motion implies that 

E"^™ as a function of the positions and velocities satisfies the identity 

2PN r 1 in 5^2PN- 



+ 



0(6) 



(B7) 



from which we find the law of variation of E"^^"^ under the complete 2.5PN dynamics of the 
binary, namely 

(B8) 



dt ^ & ^ dv\ 

This equation makes it clear how the damping acceleration term D\ drives the variation of 
the energy (see Eq. (15.3) in Ref. ^|). Since the right side of ()B8|) is a small post-Newtonian 
term it can be evaluated by inserting in place of E"^^^ its Newtonian value given by the first 
two terms in ()B6|) . Thus 

dt — ' c 

Using the expression ()B2d|) of the damping term we arrive after a short calculation at the 
balance-like equation 



5^4"^it;iDi + 0(6) . 



(B9) 



dE 
~dE 



2.5PN 



-£^ + 0(6) 



(BIO) 



which involves in the left side the 2PN energy (B6) augmented by a pure 2.5PN term, namely 



E 



2.5PN 



8 G^m^z/^ 



[nv IV 



(Bll) 
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and in the right side the expression 



C^ = ^9!l!plll2v'^n{nvf] . (B12a) 



15c5 

It is simple to re-write the latter expression in the well-known form 

^'' = -^Q^fQ'if+0i7). (B12b) 

As this is the standard (Newtonian) quadrupole formula we conclude that Eq. ()B10|) proves 
the energy balance between the loss of orbital energy .E^.spn ^^iq binary and the Newtonian 
energy flux carried out by the gravitational waves. Note that it is important that the energy 
which enters the left side of the balance equation is given as an instantaneous functional of 
the two world- lines of the binary (see Damour ll| for a discussion). 
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